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Abstract 

Using the SANC system we study the one-loop electroweak standard model pre- 
diction, including virtual and real photon emissions, for the decays of on-shell vector 
and scalar bosons B — > ff(j), where B is a vector boson, Z or W, or a Standard 
Model Higgs. The complete one-loop corrections and exact photon emission matrix 
element are taken into account. For the phase-space integration, the Monte Carlo 
technique is used. For Z decay the QED part of the calculation is first cross-checked 
with the exact one- loop QED prediction of KORALZ. For Higgs boson and W decays, a 
comparison is made with the approximate QED calculation of PHOTOS Monte Carlo. 
This provides a useful element for the evaluation of the theoretical uncertainty of 
PHOTOS, very interesting for its application in ongoing LEP2 and future LC and 
LHC phenomenology. 
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1 Introduction 



The SANC project of Ref. M has several purposes. The intermediate goal is to summarize 
and consolidate the effort of the last three decades in calculating Standard Model radiative 
corrections for LEP, in a well organized calculational environment for future reference. 
However, it is aimed not only at training young researchers and students, but at some 
remaining calculational projects for LEP as well. 

SANC provides an Internet-oriented, graphic interface platform, which is meant to serve 
as a starting point for longer term research efforts in the area of the higher order (multi- 
loop) calculations within the SM and beyond, for the experiments at future high energy 
colliders. An important lesson from the LEP experiments || is that the desirable way 
of providing theoretical predictions is in the form of Monte Carlo event generators. This 
aspect has been taken into account in the development of SANC from an early stage of its 
development. 

The currently available version of SANC can construct one-loop spin amplitudes for the 
decays of the gauge bosons W and Z and Higgs boson H. All of the Born and one-loop- 
corrected spin amplitudes are generated interactively from scratchP] by SANC with the help 
of the algebraic package Form3 J4j in the form of Fortran77 source codes, and are then 
used in the MC generation/intgration part of the package. For the moment, SANC features 
single real photon emission, in the calculations of the total rate and decay spectra of the 
B — > ff(j) process. The complete spin polarization density matrix of the decaying boson 
is taken into account as well. 

The integration, with the Monte Carlo method, over the three (two)-body final state 
is done without any approximation (in particular small mass approximation is not used). 
The program provides MC events with constant weight (unweighted events). The whole 
system is, therefore, fairly self-contained and complete^. 

It is of the utmost importance for such a new system to precisely reproduce known 
results. Concerning virtual corrections, comparisons were performed with the calculations 
of refs. ||,|6|]. The precise comparisons with the FeynArt package |7|] were made in ||[8,10 
as well. 

Our paper is organized as follows: in the next section we describe in detail the set of 
observables chosen for tests and the input parameter set-ups for comparisons SANC with 
KORALZ [[TJ and SANC PHOTOS @. Section 3 is devoted to a discussion of SANC technical 
reliability in the domain of QED bremsstrahlung. To this end, comparisons with KORALZ 
and leading-log distributions from PHOTOS are collected. Section 4 describes comparisons 
of SANC with PHOTOS, focusing on the non-leading contributions in W and H decays and 
on the relevant discussion of PHOTOS physical uncertainties. The summary in section 5 
closes the paper. 



1 The form2 codes related to the review of ref. || were exploited at the early stage of the SANC 
development. 

2 That is why it may be valuable for teaching and training. 
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2 Initialization set-ups for SANC , KORALZ and PHOTOS 
runs 

In the following sections we compare predictions from the programs SANC, KORALZ and 
PHOTOS. It is essential that the initialization be identical in all cases and close to the 
physical reality; in particular the following options are set in the three programs: 



In SANC we switch off the EW part of the radiative corrections. The soft photon 
limit is kept at 0.005 of the decaying particle mass. 

• In KORALZ we switch to the 0(a) mode of operation (no exponentiation) with final 
state bremsstrahlung only. We set the CMS energy equal to the Z mass. The s- 
channel 7 exchange is switched off and the soft photon limit is kept at 0.01 of the 
"beam energy" . 

• In PHOTOS we switch off the double bremsstrahlung corrections but we keep inter- 
ference effects. The soft/hard photon limit is kept at 0.005 of the decaying particle 
mass. For the generation of the Born-level two-body decays, we use the Monte Carlo 
generation from SANC. 

To visualize the differences (or the agreement) between the calculations, we choose 
a certain class of (pseudo-)observables, more precisely the one- dimensional distributions, 



which are quite similar to the one used in the first tests of PHOTOS reported in ||13|| . To 
visualize the usually small differences, we plot ratios of the predictions from two programs 
rather than the distributions themselves. 

List of observables: 

• -A- Photon energy in the decaying particle rest frame: this observable is sensitive 
mainly to the leading-log (i.e. collinear) non-infrared (i.e. not soft) component of 
the distributions. 

• -B- Energy of the final state charged particle: as in the previous case this observable 
is sensitive mainly to the leading-log (i.e. collinear) non-infrared (i.e. not soft) 
component of the distributions. 

• -C- Angle of the photon with respect to one of the charged final state particles: this 
observable is sensitive mainly to the non-collinear (i.e. non-leading-log) but soft 
(i.e. infrared) component of the distributions. 

• -D- Acollinearity angle of the final state charged particles: this observable is sensitive 
mainly to the non-collinear (i.e. non-leading-log) and non-soft (i.e. non-infrared) 
component of the distributions. 
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3 Technical tests of SANC 



Let us first cross check SANC predictions in the case of Z — > decay and for the real 

photon emission with KORALZ. The corresponding part of the KORALZ code, essentially 
an improved emulation of the MC program MUSTRAL [[TJ[], was thorougly tested over two 
decades and is an excellent candidate for the benchmark test for SANC matrix-element 
and phase-space generation. It can be noted that the phase-space parametrization of the 
single photon emission in SANC, although developed independently, is essentially the same 

as in KORALB [H, 101 ED. 

As we can see from figs. 1 and 2 the level of the agreement for all distributions of 
the types A, B, C and D is satisfactory, at the level of better than 0.5% (or at the level 
of statistical error). We have checked that the 0(a) FSR correction to the total decay 
rate calculated as a difference of SANC results, with FSR on and off, agrees well with the 
standard factor 1 + f f = 1.001743. In fact we get 1.001733 ± 2.8 x 10~ 5 . The agreement 
is thus better than 10~ 4 . We can thus conclude that SANC is "commissioned" for the Z 
decay. 
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Figure 1: Comparisons (ratios) of the SANC and KORALZ predictions for the Z decay. 
Observables A and B: ratios of the photon energy (left-hand side) and muon energy (right- 
hand side) distributions from the two programs. The dominant contribution is of leading- 
log (collinear) nature. 



In the case of the W and Higgs decay we rely on the comparison with PHOTOS. Due to 
the incomplete one-loop QED in PHOTOS we cannot, in priciple, expect the agreement for 



the decay distributions A, B to be better than the order of 



lnm|/mj 



7%. We see in fig. 



3 for Higgs boson decay and in fig. 5 for W decay, that this is indeed the case. In fact, 
the agreement is much better in the case of the Higgs boson decay. In figs. 4 and 6 we see 
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Figure 2: Comparisons (ratios) of the SANC and KORALZ predictions for the Z decay. 
Observables C and D: ratios of the photon angle with respect to [f (left-hand side) and 
fj.~H + acollinearity (right-hand side) distributions from the two programs. The dominant 
contribution is of infrared non-leading-log nature for the left-hand side plot, and non- 
infrared non-leading-log nature for the right-hand side one. 



the angular distributions of the photon with respect to the charged fermion (observable 
C). Here, in principle, the agreement should not be worse than lnmB 1 / Emin — 20% for the 

photon emitted in directions far from the charged particle and lnm 2 / m 2 \ nmB /E min — 1-4% 
for directions close to the charged particle. This is indeed the case. 



4 Physical uncertainties of PHOTOS in W and H decay 

In the case of testing H and W decays, the comparisons of PHOTOS with the "matrix- 
element" type calculations were never accounted for in a well documented way. The 
following discussion will give us a chance systematically evaluate uncertainties of PHOTOS 
for these two decays. 

In the case of the Linear Collider studies and many scenarios for new physics, the 
cross section for the production of the Higgs boson can be quite large. In the case of 
sophisticated observables, the reconstruction of the reference frames of the intermediate 
states requires good control of the bremsstrahlung corrections. Our comparison shows 
that, in the case of the Higgs boson decay, the approximation used in PHOTOS works much 
better than could be expected from its design principles. In fact it provides results difficult 
to distinguish from the "matrix-element" ones in the case of all observables A to D (figs. 
3 and 4). 
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Figure 3: Comparisons (ratios) of the SANC and PHOTOS predictions for the H decay. 
Observables A and B: ratios of the photon energy (left-hand side) and muon energy (right- 
hand side) distributions from the two programs. The dominant contribution is of leading- 
log (collinear) nature. 



At LEP2, the production and decay of W pairs is now being combined for all four 
LEP2 experiments, and uncertainties due to bremsstrahlung in W decay are important. 
PHOTOS is part of one of the main programs (YFSWW3) used in the LEP2 analysis of the 
VF-pair data. Our tests will provide an estimate of the size of the uncertainties, due to use 
of PHOTOS. If they turn out to be sizeable, the improvements will need to be implemented 
either into the PHOTOS algorithm or by other means. 

In figs. 5 and 6, we see that up to leading order, the distributions agree with the 
"matrix element" results provided by SANC. There is, however, no exceptionally good 
agreement observed in the previous case of H decay. 

This point requires further study, before planning improvements of PHOTOS. In par- 
ticular, a comparison of PHOTOS with another available matrix-element calculation of 
t — > eiyu(j), e.g. from TAUOLA should also be repeated to check if a pattern of 

discrepancy similar to that for W decay is observed. The W channel is of importance for 
some studies |01 of LHC Higgs discovery potential as well, see also 



5 Summary 

We have successfully tested SANC versus KORALZ in case of the Z decay. In this way we 
have verified the technical correctness of SANC, e.g. its phase-space generation and QED 
correction amplitudes. 
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Figure 4: Comparisons (ratios) of the SANC and PHOTOS predictions for the H decay. 
Observables C and D: ratios of the photon angle with respect to [f (left-hand side) and 
fj.~H + acollinearity (right-hand side) distributions from the two programs. The dominant 
contribution is of infrared non-leading-log nature for the left-hand side plot, and non- 
infrared non-leading-log nature for the right-hand side one. 



In case of W and H decays, we have checked that SANC agrees with the leading order 
QED calculation provided by the PHOTOS Monte Carlo. These comparisons allow us 
to evaluate the size of missing non-leading terms in PHOTOS. In the case of H decays, 
we have found that PHOTOS results are exceptionally good - differences with SANC are 
indistinguishable from zero and below 1% everywhere. In the case of W decay, PHOTOS 
predictions are within 7% for the end parts of the spectra affected by the leading-log 
corrections and within 20% for the angular part of the distributions, where the infrared- 
induced logarithm dominates over non-infrared non-leading terms only. The differences 
are up to 40% in the phase-space regions where only non-leading corrections contribute 
to the matrix element. We can conclude that results from all three programs for all four 
observables agree within the expected precision ranges. 
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Figure 5: Comparisons (ratios) of the SANC and PHOTOS predictions for the W decay. 
Observables A and B: ratios of the photon energy (left-hand side) and muon energy (right- 
hand side) distributions from the two programs. The dominant contribution is of leading- 
log (collinear) nature. 
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